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ABSTRACT 

We report the first simulations of the dynamics of quantum logic operations with a large 
number of qubits (up to 1000). A nuclear spin chain in which selective excitations of 
spins is provided by the gradient of the external magnetic field is considered. The spins 
interact with their nearest neighbors. We simulate the quantum CONTROL-NOT (CN) 
gate implementation for remote qubits which provides the long-distance entanglement. Our 
approach can be applied to any implementation of quantum logic gates involving a large 
number of qubits. 
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1. Introduction 



The field of quantum computation has achieved three important milestones: the first 
quantum algorithm the first error correction code 0, and the first experimental im- 
plementation of quantum logic |§. The next promising important step is implementation 
of quantum logic in solid-state systems with large number of qubits, say 1000 qubits. It is 
not clear which system will be the most feasible for quantum computation: nuclear spins 
[|]-||, electron spins 0j-@, quantum dots [|K|, or Josephson junctions [|ll|]-[|l^]. For all 



of these implementations, the design of a quantum computer requires simulations of the 
quantum computation dynamics on a conventional digital computer to test the quantum 
computer experimental devices. 

To simulate a general quantum computation involving N qubits, one must solve time- 
dependent equations involving all 2^ states of a quantum computer. Generally, this prob- 
lem cannot be solved on the digital computer. However, it is possible to simulate quantum 
logic operations which involve a limited number of states. These simulations can give in- 
sight into the dynamical properties of a quantum computer. Simulations of experimental 
implementations of quantum logic operations can explore the advantages and disadvan- 
tages of experimental devices long before they are built. In this paper, we report the first 
simulation of quantum logic for a large number of qubits (up to 1000). In Sec. 2, we 
describe the nuclear spin quantum computer which we simulate. In Sec. 3, we consider 
the Hamiltonian of the nuclear spin chain and the equations of motion for the amplitudes 
of the quantum states. In Sec. 4, we discuss resonant and non-resonant transitions in 
the spin chain under the action of radio- frequency (rf) pulses. In Sec. 5, we describe a 
quantum CN gate which entangles the two qubits at opposite ends of the spin chain. In 
Sec. 6, we give analytical analysis of the CN gate. In Sec. 7, the results of our simulations 
are presented. In the Conclusion we summarize our results. 

2. Nuclear spin quantum computer 



We consider a chain of identical nuclear spins placed in a high external magnetic field, 
B (Fig. 1). We suppose that these spins are initially polarized along the direction of the 
external field (z-direction). The NMR frequency is fo = (7/27r)I?o, where 7 is the nuclear 
gyromagnetic ratio. For example, for a proton in the field B = 10T, one has the NMR 
frequency f « 430MHz. 

Next, we suppose that the external magnetic field is slightly non-uniform, Bq = B$(z). 
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Suppose that the frequency difference of two neighboring spins is, Af « 10kHz. Thus, if 
the frequency of the edge spin is 430MHz, the frequency of the other edge spin is ~ 440MHz. 
Then, the value of Bo increases by ABq = 0.23T along the spin chain. Taking the distance 
between the neighboring spins, a m 2A, we obtain the value for the gradient of the magnetic 
field, \dBo/dz\ »s 0.23/1000CICOS6 1 , where 9 is the angle between the direction of the chain 
and the 2-axis (Fig. 1). Below we will take cos# = 1/VS. Thus, the gradient of the 
magnetic field is \dB /dz\ k2x 10 6 T/m. 

Next, we discuss the interaction between spins. In a large external magnetic field, B Q , 
the stationary states of the chain can be described as a combination of individual states of 
nuclear spins, for example, 

|00...00), |00...01), 

and so on, where the state |0) corresponds to the direction of a nuclear spin along the 
direction of the magnetic field and the state 1 1) to the opposite direction. The magnetic 
dipole field on nucleus j in any stationary state is much less than the external field. So, 
only the ^-component of the dipole field, B dz = B dz (j), affects the energy spectrum, 

BdzW = 2^ 3 Vkz, {ktJ), I 1 ) 

k=0 T kj 

where fi^z is the ^-component of the nuclear magnetic moment, is the distance between 
the nuclei k and j. To suppress the dipole interaction, one should choose the angle 9 m 
54.7°, for which cos9 = l/y/3. Then, for any stationary state, the ^-component of the 
dipole field disappears. 

We assume that main interaction between the nuclear spins (when the dipole interaction 
is suppressed) is an Ising type of interaction mediated by chemical bonds. This situation 
is observed in liquids where the dipole-dipole interaction is suppressed by the rotational 
motion of the molecules. Nuclear spins in liquids were used for quantum computations 
involving a small number of qubits [14]| . 



3. The Hamiltonian and equations of motion 



The Hamiltonian for the chain of spins considered can be written in the form [15 

N-l N-2 

h = -y, - 2 j e nn+i + v, (2) 

k=0 k=0 
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where uj^ is the Larmor frequency of the fc-th spin (neglecting interactions between spins), 
Dk = r yB (zk), J is the constant of the Ising interaction, 7| is the operator of the z- 
component of spin 1/2, the operator, V, describes the interaction with pulses of the rf 
field, Zk is the ^-coordinate of the fc-th spin, and we set % = 1. 

Below, we assume that the characteristic values of the parameters in the Hamiltonian 
(2) are, 

uj k /27T = f + kAf, f w 430MHz, Af « 10kHz, J/2n « 100^, iV < 1000. (3) 
The operator V for the n-th rf pulse can be written |0 



n(n) AT-1 
z fc=0 L 



4" exp(-iu; (n) t) + 4+ exp(wj (n) t) 



(4) 



where is the Rabi frequency of the n-th pulse, I k = 7|?±i/f , and u/ n > is the frequency 
of the n-th pulse. We choose the value of Q ~ 0.1J. (We assume that the rf field is circularly 
polarized in the xy plane.) 

In the interaction representation, the wave function, of the spin chain can be written 

as, 

^ = Y, C P\P) exp(-z^), 
p 

where E p is the energy of the state \p). Substituting the expression for the wave function 
^ into the Schrodinger equation, we obtain the equation of motion for the amplitude C p , 

2 N -1 

iC p = ]T C ™V$ exp[i(E p - E m )t + ir pm ^ n H}, (5) 

m=0 

where r pm = =)=1 for E p > E m and E p < E m , respectively, = —VL^/2 for the states 
\p) and \m) which are connected by a single-spin transition, and = for all other states. 

4. Resonant and non-resonant transitions 

As the number of spins, N, increases, the number of states increases exponentially but 
the number of resonant frequencies in our system is 3iV — 2 because only single-spin tran- 
sitions are allowed by the operator V in (4). The resonant frequencies of our spin chain 

are, 

cu k ± J, (k — 0, or k = N - 1), (6) 
u k , u k ±2J,(l < k < N — 2). 
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For edge spins with k — and k = N — l, the upper and lower signs correspond to the states 
|0) or |1) of the only neighboring spin. For inner spins with 1 < k < N — 2 the frequency 
Uk corresponds to having nearest neighbors whose spins are in opposite directions to each 
other. The "+" sign corresponds to having the nearest neighbors in their ground state. 
The "-" sign corresponds to having the nearest neighbors in their excited state. 
Now we consider any basic stationary state, 

\q N -iqN-2--qiqo), (7) 

where the subscript indicates the position of the spin in a chain, and — 0,1. If one applies 
to the spin chain a resonant rf pulse of a frequency, uj, from (6) one has two possibilities: 

1) . The frequency of the pulse uj is the resonant frequency of the state (7). 

2) . The frequency u differs from the closest resonant frequency of the state (7) by the 
value 2J or 4J. 

In the first case, one has a resonant transition. For the second case, one has a non- 
resonant transition. If J <C 27r A/ we can neglect all other non-resonant transitions for the 
state (7). Below, we will write a rigorous condition for Af which is required in order to 
neglect all other non-resonant transitions. 

Thus, considering the transformation of any basic state under the action of an rf pulse 
with a frequency, uj, from (6), we should take into consideration only one transition. This 
transition will be either a resonant one or a non-resonant one with the frequency difference 
2J or 4J. 

This allows us to simplify equations (5) to the set of two coupled equations, 

iC p = exp[t(E p - E m - uj {n) )t]C m , (8) 

iC m = exp[i(E m -E p + J n) )t]C p , 

where E p > E m , \p) and \m) are any two stationary states which are connected by a 
single-spin transition and whose energies differ by a/") or uj^ ± 2 J or uj^ ± 4J. 

The solution of equations (8) for the case when the system is initially in a stationary 
state \m), can be written, 

C m (t + t) = [cos(fi e r/2) + i(A/Sl e ) sin(fi e r/2)] x exp(-irA/2), (9) 

C p (t + r) = i(Q/Q e ) sin(fi e r/2) x exp(it A + it A/2), 
C m (t ) = 1, C p (t ) = 0. 

In Eqs (9) we omitted the upper index "n" which indicates the number of the rf pulse, t is 
the time of the beginning of the pulse, r is its duration, A = E p — E, m — uj, il e = (f2 2 + A 2 ) 1 / 2 
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is the NMR frequency in the rotating frame. If the system is initially in the upper state, 
\p)i (C m (t ) = 0, C p (t ) = 1), one can obtain the solution of Eqs (8) by changing the sign 
at A and setting: m — > p and p — > m in (9). 

For the resonant transition (A = 0) the expressions (9) transform into the well-known 
equations for the Rabi transitions, 

C m {t + r) = cos(fir/2), C p {t + r) = i sin(fir/2). (10) 

In particular, for VLt = 7T (a 7r-pulse), Eqs (10) describe the complete transition from the 
state \m) to the state \p). 

For non-resonant transitions, expressions (9) include two characteristic parameters: 
Q/Q e and sin(fi e r/2). If either of these two parameters is zero, the probability of a non- 
resonant transition disappears. The second parameter is equal zero when Q e r = Ink 
(k = 1, 2, ..), where k is the number of revolutions of a non-resonant (average) spin about 
the effective field in the rotating frame. This is the basis of the "2-7r/c" -method for elimi- 



nation non-resonant transitions. (See ||15|| , Chapter 22, and |Tj| 



5. A Control-Not gate involving remote qubits and their long-distance en- 
tanglement 

A pure Control-Not {CN ab ) gate is a unitary operator which transforms the basic state, 

\qN-i--Qa <7b—<7o), 

into the state, 

\0N-i-Qa ?b—?o), 

where g& = 1 — if q a — 1; and % = Qb if Q a = 0. The a-th and 6-th qubits are called the 
control and the target qubits of the CN ab gate. A modified CN gate performs the same 
transformation accompanied by phase shifts which are different for different basic states 



[15]. It is well-known that the CN gate can produce an entangled state of two qubits, which 
can not be represented as a product of the individual wave functions. 

We shall consider an implementation of the CN gate in the Ising spin chain with the 
left end spin as the control qubit and the right end spin as the target qubit, i.e. CN^-ip 
for a spin chain of 200 and a spin chain of 1000 qubits. Using this gate we will create 
entanglement between the end qubits in the spin chain. We start with the ground state. 
Then we apply a 7r/2-pulse with frequency u^-i to produce a superpositional state of the 
(N - l)-th (left) qubit, 

* = |0...0)+i|1...0). (11) 



6 



(Here and below the normalization factor l/y/2 is omitted.) To implement a modified 
CN N _i j0 gate we apply to the spin chain L = 397 7r-pulses if N = 200, and L = 1997 
pulses if N = 1000. The first 7r-pulse has the frequency uj = ujn-2- For the second 7r-pulse 
uj = ujn-z- For the third 7r-pulse, uj = ujn-2 — 2J, etc. 

6. An analytic solution 

An analytical expression for the wave function, after the action of ir/2- and L 7r-pulses, 
can be easily derived if for all 7r-pulses Q e r/2 = 2nk, with the same value k: 

* = C o |00..0) + Ci|10...1), (12) 

C = (-l) feL exp(-mLv / 4P - 1/2), d = -1, (13) 

For k >> 1, we get the same solution for odd and even k: Co ~ 1. This result is easy 
to understand. For a 7r-pulse, the Rabi frequency is f2 = \A\/\/4k 2 — 1. Large values 
of k correspond to small values of the parameter f2/|A|. If fi/|A| approaches zero, the 
non-resonant pulse cannot change the quantum state. 

For a small value of k, the non-resonant pulse can change the phase of a state. For 
example, for k = 1 we have, 

C = exp[i7rL(l - VS/2)]. 

After the first 7r-pulse, the phase shift is approximately 24°, but it grows as the number of 
7r-pulses, L, increases. 

Now, we shall mention an important point. If we consider the probability of non- 
resonant transition, the small parameter of the problem is: 

e = (fi/fi e ) 2 sin 2 (fi e r/2). 

(It follows from (9) that the expression for |C m | 2 can be written in the form: |C m | 2 = 1 — e, 
and \C P \ 2 = e.) If we take into consideration the change of the phase of non-resonant state, 
the small parameter of the problem is fi/||A| ~ VL/Vt e . 

Next, we will discuss the probability of non- resonant transitions using perturbation 
theory. The analytical expression for probabilities |C | 2 and |Ci| 2 can be easily found in 
the first non-vanishing approximation of perturbation theory: 

|C | 2 = l-Le, |d! 2 = 1. (14) 



7 



The decrease of the probability | Co | 2 is caused by the generation of unwanted states. One 
can see that the deviation from the value | C 1 2 = 1 accumulates when the number of ir- 
pulses, L, increases. It means that the small parameter of the problem is Le rather than 
e. Consider first those non-resonant transitions which we ignore in this paper. For the 
number of qubits, iV = 1000, using the characteristic parameters from (3) and L xs 2000, 
VI ks 0.1 J, |A| = 2nAf, we obtain, Le < 10~ 3 . Thus, as was already pointed out, one can 
neglect non-resonant transitions whose frequency differences are of the order of Af. 

Now we consider non-resonant transitions which are included in our simulations. Set- 
ting |A| = 2J, we obtain: Le < 2.5. Thus, the deviation from the 2nk condition can 
produce large distortions from the desired wave function. To study these distortions we 
have to use computer simulations. 

7. Computer simulations 

We have developed a numerical code which allowed us to study the dynamics of quan- 
tum states with the probabilities no less than 10~ 6 for a spin chain with up to 1000 qubits. 
(Because we omitted in (11) the normalization factor l/y/2 all probabilities in this section 
including in Figs 2-14 are doubled.) The sum of the probabilities of all these states was 
equal 2 — O(10~ 6 ) (the normalization condition). 

Next, in Figs 2-14, we present the results of computer simulations with N = 200 and 
with N = 1000 qubits. Fig. 2, shows the probability of the excited unwanted states after 
implementation of the CNigg^ gate, for N = 200 and Q = 0.14. On the horizontal axis the 
unwanted states are shown in the order of their generation. A total of 7385 unwanted states 
were generated which had the probability, P > 10~ 6 . (In all Figs 2-14 only the states with 
P > 10 -6 are taken into account.) The probability distribution of unwanted states clearly 
contains two "bands" . One group of these states has the probability, P ~ 10~ 6 (the bold 
"line" near the horizontal axis). The second group of states has the probability P ~ 10~ 3 
(the upper "curve" in Fig. 2). Fig. 3, shows an enlargement of the upper "band" of the Fig. 
2. One can see some sub-structure of this "band". Fig. 4, shows the sub-structure in the 
lower "band" of Fig. 2. Figs 5 (a,b) show the sub-structure of the upper and lower "bands" 
shown in Fig. 4. One can see some hierarchy in the structure of the distribution function of 
generated unwanted states. Figs 6 (b-j ) show the typical structure of unwanted states of the 
spin chain. Fig. 6a, shows the ground state of the spin chain (all qubits are in their ground 
state). The value of P in Fig. 6, indicates the probability of the states. All states in Figs 
6 (b-j ) belong to the upper "band" shown in Fig. 2, i.e. they have a probability, P ~ 10~ 3 . 
It is interesting to note that the group of unwanted states with high probabilities contains 
the high energy states of the spin chain (many-spin excitations). Typical unwanted states 
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of the lower "band" in Fig. 2 (with P ~ 10~ 6 ) are shown in Figs 7(a-j) and in Figs 8(a-j). 
It is important to note that typical unwanted states for both groups (Figs 6-8) contain 
highly correlated spin excitations. Fig. 9 shows the total number of unwanted states (with 
probability P > 1CT 6 ) and the probability of the ground state, |C | 2 , as a function of the 
Rabi frequency, Q. The maximum value of |Co| 2 and the minimal total number of unwanted 
states correspond to the values of Q which satisfy the 27r/c-condition for 396 of the total 
number of 7r-pulses, 397. (The third 7r-pulse does not satisfy a 27r/c-condition.) 

Next, we have studied the generation of unwanted states for the case when only a group 
of pulses had values of Rabi frequency which deviated from the 27r/c-condition. We changed 
the value of for all 7r-pulses from k\ = 10 to k 2 = (10 + Ak). Fig. 10, shows the number of 
unwanted states and the probability of the ground state, | Co | 2 , as a function of the number 
of distorted pulses, Ak, for N = 1000. The value Q ~ 0.100 in Fig. 10 corresponds to the 
27rA;-condition for all pulses (except the 3-rd 7r-pulse), for distorted Ak pulses, Q = 0.101. 
Figs 11 and 12 demonstrate the same quantities (the number of unwanted generated states 
and |C | 2 ) for the case when Q is a random parameter for a group of pulses. (Both figures 
show typical realizations for definite distributions of Q.) Fig. 13, shows the dependence of 
the number of unwanted states and | Co | 2 on the location of the group of distorted pulses. 
Fig. 14, demonstrates the ground state and the examples of unwanted states generated 
due to the distortion of this group of pulses. One can see again that high-energy states of 
the spin chain, with many-spin excitations, are generated. 

Conclusion 

In this paper we presented the results of simulations of quantum Control-Not gate, CN N _ lj0 , 
between remote qubits, (N — l)-st and 0-th, and the creation of long-distance entangle- 
ment in nuclear spin quantum computer having a large number of qubits (up to 1000). A 
considered quantum computer is a one-dimensional nuclear spin chain placed in a slightly 
non-uniform magnetic field, and oriented in such direction that the dipole interaction be- 
tween spins is suppressed. So, the Ising interaction comes into play, as in the case of the 
liquid NMR. We used two essential assumptions: 

1. The nuclear spin chain is prepared initially in the ground state. 

2. The frequency difference between the neighboring spins due to the inhomogeneity of the 
external magnetic field is much large than the Ising interaction constant. 

Using these assumptions, we developed a numerical method which allowed us to simu- 
late the dynamics of quantum logic operations taking into consideration all quantum states 
with the probability no less that P = 10 -6 . For the case, when the 27r/c-condition is satisfied 
(the 7r-pulse for the resonant transition is at the same time a 27rA;-pulse for non-resonant 
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transitions), we obtained an analytic solution for the evolution of the nuclear spin chain. 
In the case of small deviations from the 27r/c-condition, the error accumulates. So, the 
perturbation theory becomes invalid even for small deviations from the 27r/c-condition. In 
this case, the numerical simulations are necessary. 

The main results of our simulations are the following: 

1. The unwanted states exhibit band structure in their probability distributions. There 
are two main "bands" in the probability distribution of unwanted states. The unwanted 
states in these "bands" have significantly different probabilities, Pi ow /Pupper ~ 1CT 3 . Each 
of these two bands have their own structure. 

2. A typical unwanted state is a state of highly correlated spin excitations. An important 
fact is that the unwanted states with relatively high probability include high energy states 
of the spin chain (many-spin excitations). 

3. The method developed allowed us to study generation of unwanted states and the 
probability of the desired states as a function of the distortion of rf pulses. This can be 
used to formulate the requirements for acceptable errors in quantum computation. 

The results of this paper can be used to design experimental implementations of quan- 
tum logic operations and to estimate (benchmark) the quality of experimental quantum 
computer devices. Our approach can be extended to any types of quantum computer. 
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Figure captions 

Fig. 1: Nuclear spin quantum computer (the ground state of nuclear spins). B is the 
permanent magnetic field; B\ is the radio-frequency field. The chain of spins makes the 
angle 9 with the direction of the field B . 

Fig. 2: Probabilities of unwanted states. The total number of qubits: iV = 200; Q = 0.14. 
The number of unwanted states with probabilities \C n \ 2 > 10~ 6 is 7385. The states are 
presented in the order of their generation. 

Fig. 3: The upper band of Fig. 2, shown in a larger scale. 

Fig. 4: The lower band of Fig. 2, shown in a larger scale. 

Fig. 5a: A sub-structure of the upper band of Fig. 4, shown in a larger scale. 

Fig. 5b: A sub-structure of the lower band of Fig. 4, shown in a larger scale. 

Fig. 6: (a) The ground state of the spin chain; (b-j) The typical unwanted states with 
probabilities ~ 10~ 3 . Horizontal axis shows the position of a qubit in the spin chain of 
N = 200 spins. The vertical axis shows the state |0) or |1) of the qubit. 

Fig. 7: Examples of "low energy" unwanted states from the lower band in Fig. 2. 

Fig. 8: Examples of "intermediate energy" unwanted states from the lower band in Fig. 2. 

Fig. 9: (a) Probability, |C | 2 , as a function of Q. The total number of qubits N = 200; 
(b) The total number of unwanted states. 

Fig. 10: Dependence of the total number of unwanted states as a function of AA; = k 2 — ki 
(ki = 10, N = 1000); (a) The total number of unwanted states; Q = 0.1 for all 7r-pulses 
except for the 7r-pulses with numbers k in the range: ki < k < k 2 for which Q = 0.101; (b) 
The probability |C | 2 , for the same parameters as in (a). 

Fig. 11: (a) The probability | Co | 2 , as a function of parameter e ; N = 1000; (b) The 
number of unwanted states as a function of parameter e ; Q = 0.1 for all pulses from 10-th 
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to (10 + A/c)-th for which £1 = 0.1 + e; e is a random parameter: e £ [—e , e ]. 

Fig. 12: (a) The probability, | Co | 2 , as a function of Ak = k<i — k\\ (b) The total number of 
unwanted states; Q = 0.1 for all 7r-pulses except for the 7r-pulses from 10-th to (10 + Ak)- 
th for which f2 = 0.1 +e; e is a random parameter in the range —0.05 < e < 0.05, N = 1000. 

Fig. 13: (a) The number of unwanted states as a function of the parameter k±; Q = 0.1 for 
all 7r-pulses except for the 7r-pulses from fci-th to (ki + 15)-th for which Q = 0.1 + e; e is a 
random parameter, e £ [—0.005, 0.005]; (b) The probability, |C | 2 , for the same parameters 
as in (a); N — 1000. 

Fig. 14: (a) The ground state of the chain; (b-j) Examples of unwanted states (N = 1000); 
Q = 0.1 for all 7r-pulses except for the 7r-pulses from 10-th to 40-th for which Q = 0.1 + e, 
where e is a random parameter in the range: —0.05 < e < 0.05 . 
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